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INTRODUCTION 


1.  Pavement  Evaluation 

This  document  reports  research  that  is  part  of  a  larger  Air  Force 
Weapons  Laboratory  project  lO  develop  a  pavement  evaluation  procedure.  The 
following  steps  are  involved  in  the  pavement  evaluation  procedure,  which  for 
the  present  uses  linear  elastic  analysis  for  determination  of  the  stresses  and 
strains  in  the  pavement  structure. 

The  first  step  is  to  measure  the  shear  modulus  for  different  layers  of 
the  pavement  structure  in-iitu,  using  a  nondestructive  vibration  testing  method. 
The  shear  modulus  thus  measured  will  be  for  very  small  strain  amplitudes,  on 
the  order  of  10“5  in/in  or  less.  Because  the  stress-strain  relations  for  paving 
materials  are  nonlinear,  the  secant  shear  modulus  for  larger  strain  amplitudes 
produced  by  an  aircraft  loading  will  be  smaller  than  the  modulus  measured  by 
the  nondestructive  testing  method . 

The  second  step  is  to  make  the  proper  reduction  in  the  measured 
shear  modulus  to  correspond  to  the  strain  produced  by  an  aircraft  loading. 

The  third  step  is  to  use  this  shear  modulus  in  a  finite-element  analysis 
for  stresses  and  strains  in  pavement  structure  under  load. 

The  fourth  and  final  step  is  to  assess  the  amount  of  damage  or  pave¬ 
ment  distress  that  will  be  produced  by  the  loading. 

This  report  presents  a  procedure  for  the  second  step,  the  reduction 

of  shear  modulus  with  increasing  strain  level. 

1 
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2.  Review  of  Previous  Report 

The  first  phase  in  developing  a  procedure  for  the  reduction  of  shear 
modulus  with  increasing  strain  level  was  to  design  and  construct  testing  equip¬ 
ment  capable  of  accurate  measurement  of  the  shear  stress-strain  relation  for 
soils,  over  a  wide  range  of  strain  levels.  The  range  was  from  about  10”^  in/in 
to  failure,  i.e. ,  after  a  given  number  of  constant-amplitude  cyclic  loads  each 
sample  was  loaded  to  failure.  Also,  a  series  of  tests  with  two  soils,  a  silty 
sand  and  a  silty  clay,  were  conducted  to  assess  the  relative  effects  of  various 
parameters,  such  as  density,  percent  saturation  anu  confining  stress,  on  the 
shear  stress-strain  relation  for  cyclic  loading.  Reference  1  reports  on  this 
phase  of  the  research,  and  the  testing  equipment,  its  capabilities,  and  the  testing 
procedure  developed  are  described  in  detail.  Examples  of  recorded  stress-strain 
curves  are  given,  and  the  methods  of  analysis  of  the  data  are  discussed.  The 
same  equipment,  procedures,  and  most  of  the  methods  of  analysis  were  used  for 
the  current  testing  program.  They  will  not  be  deset  ibed  in  detail  herein,  since 
the  information  is  available  in  reference  1.  However,  a  brief  description  of 
methods  and  procedures  is  given  in  Appendix  n.  Some  of  the  definitions 
presented  in  Section  n  are  amplified  in  reference  1. 

3.  Current  Testing  Program 

The  objective  of  the  current  testing  program  was  to  determin'  whether 

1.  Hardin,  Bobby  O. ,  Characterization  and  Use  of  Shear  Stress-Strain 

Relations  for  Airfield  Subgrade  and  Base  Course  Materials,  Technical 
Report  No.  AFWL-TR-71-60,  Air  Force  Weapons  Laboratory,  KirtJand 
AFB,  NM,  July  1971. 
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or  not  the  characterization  of  the  shear  stress-strain  relation  presenter!  in 
reference  1  is  applicable  to  a  wide  variety  of  soils.  A  total  of  80  specimens  of 
20  different  soils  were  tested  with  differing  parameters,  such  as  soil  type  , 
strain  amplitude,  density,  and  percent  saturation.  Including  Phase  I,  a  total  of 
129  tests  on  24  different  soils  have  been  conducted  and  analyzed.  Data  for  only 
six  tests  appeared  to  be  faulty  (i.  e. ,  an  error  in  conducting  the  test)  and  were 
thrown  out.  The  relationships  presented  here  arc  based  on  123  tests  of  24 
different  soils. 

A  procedure  for  reducing  the  shear  modulus  with  increasing  strain  level, 
recommended  for  use  in  the  Air  Force  pavement  evaluation  technique,  is  pre¬ 
sented  in  Section  II.  In  Section  III  values  of  the  shear  modulus  determined  by 
this  procedure  are  compared  to  the  experimental  data  in  order  to  verify  the 
procedure  and  to  show  the  probable  magnitude  of  error  in  using  the  procedure. 

A  few  of  the  tests  involved  more  general  loading  histories,  with  mixed  loading 
amplitudes  and  rest  periods  between  loads.  Results  and  discussions  of  these 
tests  are  presented  in  Section  IV. 
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SECTION  II 


A  PRACTICAL  PROCEDURE  FOR  REDUCTION  OF  SHEAR  MODULUS 
WITH  INCREASING  STRAIN  LEVEL 


1.  Objective 

The  objective  of  this  section  is  to  present  the  procedure  recommended 
for  reduction  of  shear  modulus  with  increasing  strain  level  in  a  practical  form, 
for  use  in  pavement  evaluation,  unencumbered  by  details  of  testing  or  presen¬ 


tation  of  the  supporting  data. 
2.  Definitions 


The  following  parameters  are  used  and  are  defined  with  reference  to 


figure  1: 


Maximum  Shear  Modulus  =  G  =  the  initial  tangent  modulus, 

max 


-F> 


or  secant  modulus  for  strain  amplitude  £  10  in/in  (For 


pavement  evaluation  this  quantity  is  to  be  measured  by  the 
nondestructive  vibratory  test); 


Maximum  Shear  Stress  -  T  strength  of  the  specimen  in 

max 


simple  shear,  defines  assymptote  in  figure  1  (This  quantity 


has  been  related  emperically  to  G  see  Appendix  I. 


However,  if  measured  values  of  r  are  used,  a  rate  of 

max 


.oading  corresponding  to  T  -  0. 1  will  suffice,  eventhough 
the  value  of  T  for  actual  aircraft  loading  may  be  0.01); 


Reference  Strain  y  ~  r  /G  ,  defined  by  the  intersection  of 

r  max  max 


the  initial  tangent  line  and  strength  assymptote  in  figure  1, 


and  emperically  related  to  G  in  this  section; 

max 


Normalized  Strain  y/y  ,  where  y  i.s  the  shear  strain; 

Hyperbolic  Strain  y^  a  parameter  defined  by  equation  9  {If  equation  9 
is  substituted  into  equation  8,  the  relationship  between  normalized 
shear  modulus  and  normalized  strain  is  obtained.  This  relation¬ 
ship  changes  with  the  value  of  a,  as  shown  in  figure  15,  depending 
on  the  values  of  S,  N,  T,  and  the  type  of  soil;  see  equation  19. 

An  alternate  method  of  presentation  is  used  in  this  section. 

Graphs  based  on  equation  9  are  presented  giving  the  value  of  y^ 
corresponding  to  a  value  of  y/y  ,  depending  on  the  values  of  S, 

N,  T,  and  soil  type.  With  the  value  of  y^  determined,  the 
normalized  shear  modulus  is  given  by  the  simple  hyperbolic 
equation  8.); 

Shear  Modulus  =  G  =  secant  shear  modulus  for  a  given  strain  amplitude 
and  cycle  as  shown  in  figure  1; 

Normalized  Shear  Modulus  =  G/G  ; 

max  ’ 

Number  of  Cycles  -  N  ; 

Strain  Time  T  time  in  minutes  to  reach  a  normalized  strain 
equal  to  one; 

Percent  Saturation  -  S  =  ratio  of  volume  of  water  to  volume  of  voids 
in  the  soil,  expressed  as  a  percentage,  and; 

Void  Ratio  =  e  -  ratio  of  volume  of  voids  to  volume  of  solids  in  the  soil. 

Determination  of  the  Reference  Strain 

The  reference  strain  is  given  by  the  following  equation  derived  in 
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Appendix  I. 


umax 

FiF 


0.0  -  0.25  (PI) 


(1) 


where  PI  =  plasticity  index,  Gmax  is  in  psi, 


F  = 


(2.973  -  e) ' 


(1+e) 

R  -  1100,  for  sands  with  less  than  15  percent  fines 
and  R  -  1100  -  6S,  for  cohesive  soils  with  more  than  15  percent  fines 
Let 


(2) 


F2!! 


2 


1 


(4) 


0.6  -  0.25  (PI) 


076 


Then  for  sands  with  less  than  15  percent  fines,  PI  -  0  and  R  =  1100 
giving 

Cj  =  (2.017  x  105)  F2  (5) 

For  cohesive  soils  with  more  than  15  percent  fines,  substituting  R  =  1100  -  OS 


into  equation  4 


Finally 


F2  (1100  -  OS)2 
C1  ~  0.6  -  0.25  (PI)^*6 


max 


(6) 


(7) 


With  Gmax  in  psi. 

The  vi  lue  of  C-^  for  sands  with  less  than  15  percent  fines,  as  defined  by 
equations  2  and  5,  is  given  in  figure  2.  The  value  of  void  ratio  is  all  that  is 
needed  to  determine  Cj  for  this  case.  The  value  of  for  cohesive  soils  with 
more  than  15  percent  fines,  as  defied  by  equations  2  and  (3,  is  given  in  figure  3. 


For  this  case  the  values  of  the  plasticity  index  and  percent  saturation  are  needed 


in  addition  to  the  value  of  the  void  ratio.  The  dashed  line  in  figure  3  shows  how 
to  use  the  figure.  First,  using  the  left  hand  side  of  the  graph  with  the  value  of 
plasticity  index  as  abscissa,  determine  an  ordinate  depending  on  percent  saturation. 
Extend  this  ord.uate  to  the  right  to  determine  an  abscissa  for  the  right  hand  side 
of  the  graph,  corresponding  to  the  value  of  void  ratio.  The  second  abscissa 
determines  the  value  of  Cr  The  dashed  line  represents  the  case  where  PI  =  25, 

S  =  50,  and  e  •-  1.  The  value  of  Cj  --  4.9  x  10(,psi. 


Having  determined  the  value  of  C1  from  either  figure  2  or  figure  3, 
depending  on  the  soil  type,  the  reference  strain  can  now  be  determined  from 
figure  4.  Using  the  value  of  C-j  as  abscissa,  the  value  of  the  reference  strain  is 
read  as  ordinate,  from  the  curve  corresponding  to  the  value  of  Gmax  in  psi. 

A 

The  dashed  lines  in  figure  4  are  for  Cx  =  4.  9  x  10  psi  and  Gmax  -  8000  psi, 

-4 

y  =16.5x10  in/in. 
r 

4.  Determination  of  the  Shear  Modulus 

The  normalized  shear  modulus  is  given  by 

G  =  1 _ 


max 


1  +  n 


where 


1  +  a  e 


xp 


y_ 


0.4 


(8) 


(9) 


The  value  of  a  is  defined  by  one  of  the  following  equations,  depending  on  the 


type  of  soil 
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Figure  3.  Value  of  Cj  versus  plasticity  index,  percent  saturation,  and  void  ratio  for  cohesive 
soils  with  more  than  15  percent  fines. 


REFERENCE  STRAIN  IN  10  IN/ IN 


Figure  4.  Reference  strain  versus  value  of  Cj  and  maximum 
shear  modulus. 
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[(3.85/N)  -  0. 85]  T0,025 
1.6(1  +  0. 02S)  T°’  2/N°*  ° 


for  clean  dry  sands 

for  nonplastic  soils  with 
fines  and  low  plasticity  >  (10) 

soils 


[o.2(1  ,  0. OSS,  T°- 7  V- 15  T  J 

The  values  of  y^  as  defined  by  equations  9  and  10  are  given  in  figures 
5,6,7  and  8  for  various  cases.  Figure  5  is  for  clean  dry  sands.  The  solid 
curves  give  the  relationship  between  and  y/yr  for  a  fast  rate  of  loading, 

T  =  0.01.  The  dashed  curves  are  for  a  slow  rate  of  loading,  T  =  10.  Curves 
are  shown  for  N  =  1,  2,  10,  >  100.  The  dashed  and  solid  curves  coincide  for 
10  and  >  100  cycles.  Having  determined  the  reference  strain  according  to  the 
procedure  given  in  paragraph  2,  the  normalized  strain  can  be  computed,  and  the 
hyperbolic  strain  determined  from  figure  5  for  clean  dry  sands.  Similar  graphs 
are  given  in  figures  6  and  7  for  nonplastic  soils  with  fines  and  low  plasticity 
soils.  In  figures  6  and  7  the  solid  curves  are  for  60  percent  saturation  and  the 
dashed  curves  are  for  100  percent  saturation.  Figure  6  is  for  a  fast  rate  of 
loading,  T  =  0.01, and  figure  7  is  for  a  slow  rate  of  loading,  T  =  10.  Both  sets 
of  curves  in  figures  6  and  7  approach  the  hyperbolic  line  for  a  large  number  of 
cycles.  Again,  knowing  the  value  of  reference  strain,  the  value  of  the  hyperbolic 
strain  can  be  determined  from  these  figures  for  a  given  value  of  strain.  Figure 
8  is  a  similar  graph  for  high  plasticity  soils.  For  the  fast  rate  of  loading, 

T  =  9  01,  the  hyperbolic  line  defines  the  relationship  for  all  numbers  of  cycles 


and  percents  saturation.  For  the  slow  rate  of  loading,  T  =  10,  the  solid 


curves  define  the  relationship  for  60  percent  saturation,  and  the  dashed  curves 
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HYPERBOLIC  STRAIN, 


NORMALIZED  STRAIN,  J/Jr 


Fissure  5.  Hyperbolic  strain  versus  normalized  strain  for  various 
numbers  of  cycles  and  rates  of  loading. 
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HYPERBOLIC  STRAIN, 


LOW  PLASTICITY  SOILS  AND 
NONPLASTIC  SOILS  WITH  FINES 


NORMALIZED  STRAIN,  J/7r 

Figure  6.  Hyperbolic  strain  versus  normalized  strain  for  various 
numbers  of  cycles  and  percents  saturation,  slow  rate  of 
loading,  T  =  0.01,  nonplastic  soils  with  fines  and  low 
plasticity  soils. 
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HYPERBOLIC  STRAIN, 


LOW  PLASTICITY  SOILS  AND 
NONPLASTIC  SOILS  WITH  FINES 


NORMALIZED  STRAIN,  J/Jr 


Figure  7.  Hyperbolic  strain  versus  normalized  strain  for  various 
numbers  of  cycles  and  percents  saturation,  slow  rate 
of  loading,  T  =  10,  nonplastic  soils  with  fines  and  low 
plasticity  soils. 
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HYPERBOLIC  STRAIN, 


HIGH  PLASTICITY  SOILS 
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4 


3 
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0  1  2  3  4  5 

NORMALIZED  STRAIN.,  J/7r 


Figure  8.  Hyperbolic  strain  versus  normalized  strain  for  numbers 
of  cycles,  percents  saturation,  and  rates  of  loading, 
high  plasticity  soils. 
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are  for  100  percent  saturation.  With  the  value  of  hyperbolic  strain  determined 
from  one  of  figures  5,  6,  7,  or  8,  depending  on  soil  type,  the  normalized  shear 
modulus  can  be  determined  from  figure  9. 

5.  Examples  of  Use  of  the  Procedure 

For  pavement  evaluation  the  type  of  soil  (including  PI  and  particle  size) 
and  values  of  e  and  S  would  have  to  be  estimated  from  available  knowledge  of  the 
subgrade  or  from  a  core  sample;  the  value  of  G^  would  be  measured  by 
nondestructive  vibratory  testing;  the  value  of  N  would  come  from  traffic  records; 
the  value  of  y  would  be  determined  by  the  finite -element  analysis,  and;  the  value 
of  T  would  depend  on  the  speed  of  the  aircraft. 

Example  number  1: 

Given:  A  clean  dry  sand, 

Gmax  =  18380  Psi’ 

e  =  0.62, 

N  -  1 , 

y  =  18. 6  x  10~4  in/in, 

T  =  81  (a  very  slow  rate  of  loading). 

Find:  The  shear  modulus,  G; 

From  figure  2,  using  e  =  0.62;  --  2.36  x  107  psi. 

From  figure  4,  using  C-,  -  2.36  x  107  and  G  =  18380  psi; 

max 

yp  =  7.5  x  10-4  in/in. 

—4  -4 

Calculate,  using  y  =  18. 6  x  10  and  yr  =  7. 5  x  10  ; 

y/yr  =  (18.6  x  10-4)  /  (7.5  x  10"4)  =2.48. 
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From  figure  5,  using  y/yp  =  2.48,  N  =  1,  and  T  =  81; 
yh  ^  4,4 

From  figure  9.  using  yh  =  4.4;  G/G  =  0. 185 

“  max 

Calculate,  using  G/G  -  0. 185  and  G  18380  psi; 

max  max  1 

G  -  (0. 185)  (18380)  3400  psi 

The  same  parameters  given  in  this  example  were  measured  for  test  no.  101  of 
the  WES  Sand  (see  figure  16).  The  measured  shear  modulus,  G  3320  psi,  in 
this  case,  was  less  than  3  percent  differenc  ■*  between  calculated  and  measured 


values. 


Example  number  2: 

Given:  A  low  plasticity  soil, 


PI  =  6, 


Percent  passing  no.  200  sieve  =  96, 


G  =  12680  psi, 
max 

e  =  0;  67, 


S  =  73, 


N  =  10, 

y  =  12. 9  x  10~4  in/in, 


T  =  0.38  (a  medium  fast  rate  of  loading). 

Find:  The  shear  modulus ,  G; 

From  figure  3,  using  PI  =  6,  S  =  73,  and  e  =  0.67; 
C ^  =  8. 0  x  106  psi. 


From  figure  4,  using  C..  =  8.0  x  10  psi  and  G  =  12680  psi; 

X  m  aX 


yr  -  15  x  10“4  in/in. 

Calculate,  using  y  12.9  x  10-4  and  yr  15  x  10~4; 

y/yf  -  (12.9  x  lo"4)  /  (15  x  10~j  =  0.860. 

Since  T  =  0.38,  use  the  average  between  figures  6  and  7, 

also  using  y/yr  -  0.860,  N  =  10,  and  S  =  73; 

=  from  figure  6)  +  (1.44  from  figure  7)]  /2  1.22. 

From  figure  9,  using  y,  1.22;  G/C  0.450 

n  max 

Calculate,  using  G/G  0.450  and  G  12680  psi; 

max  max 

G  =  (0.450)  (12680)  5700  psi 

The  same  parameters  given  in  this  example  were  measured  for  test  no.  27  of  the 
Vicksburg  Loess.  The  measured  shear  modulus.  G  =  5580  psi,  in  this  case, 
was  less  than  3  percent  difference  between  calculated  and  measured  values.  Not 
all  cases  would  compare  this  well,  but  these  were  not  chosen  specifically  to 
show  a  good  comparison. 
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COMPARISON  OF  VALUES  FROM  THE  PRACTICAL  PROCEDURE 
TO  EXPERIMENTAL  DATA 


1.  Soils  Tested 

Classification  data  for  each  of  the  soils  tested  and  various  test  parameters 
for  each  of  the  tests  are  given  in  tables  1,2,3,  and  4.  Ir  these  tables  the  soils 
are  grouped  into  four  different  categories.  Table  1  gives  data  for  clean  sands 
and  other  sands  with  high  permeability.  Table  2  gives  data  for  nonplastic  silty 
sands.  Table  3  gives  data  for  low  plasticity  soils  and  table  4  for  high  plasticity 
soils.  The  particle  size  distribution  curves  for  these  soils  are  given  in  figures 
10,  11,  and  12.  Those  in  figure  10  are  for  nonplastic  soils,  corresponding  to 


Table  1.  Data  for  Simple  Shear  Tests  of  Clean  Sands  and  Other 
Sands  with  High  Permeability. 
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Table  2.  Data  for  Simple  Shear  Tests  of  Nonplastic  Silty  Sands. 
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Table  2.  Data  for  Simple  Shear  Tests  of  Nonplastic  Silty  Sands  (eont') 
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Table  3.  Data  for  Simple  Shear  Tests  of  Low  Plasticity  Soils. 
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Table  3.  Data  for  Simple  Shear  Tests  of  Low  Plasticity  Soils  (cont') 
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Table  4.  Data  for  Simple  Shear  Tests  of  High  Plasticity  Soils. 
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Figure  10.  Particle  size  distribution  curves  for  nonplastic  soils. 
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Figure  12.  Particle  size  distribution  curves  for  high  plasticity  soils. 
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Figure  13.  Comparison  of  the  variation  of  experimental  and  calculated 
values  of  R  with  percent  saturation,  Air  Force  Silty  Sand, 
Air  Force  Silty  Clay,  Vicksburg  Loess,  Vanceburg,  Allen, 
Kentucky  55,  Longhorn,  and  West  Virginia  Shale. 
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experimental  values.  The  weakest  comparison  is  with  two  of  the  high  plastieity 
soils,  the  Cheeks  and  Nevada  Clays.  Here  the  measured  and  calculated  values 
differ  by  a  factor  of  2.  However,  the  Louisiana  and  San  Francisco  Clays  are 
also  high  plasticity  soils,  and  they  more  nearly  conform  to  the  pattern  of  the 
other  cohesive  soils  than  to  the  Cheeks  and  Nevada  Clays.  Hence,  the  use  of 
equation  3  for  high  plasticity  soils  is  rather  uncertain  but  is  the  best  available  at 
present. 

The  measured  values  of  R  for  sands  with  less  than  15  percent  fines  are 

shown  in  table  5.  The  average  value  from  10  tests  of  dry  WES  sand  is  1102  with 

about  10  percent  scatter.  The  average  value  for  the  St.  John's  Sand  with  tests 

at  both  88  and  100  percents  saturation  is  1224,  with  the  value  for  test  no.  70 

being  excessively  high.  The  average  for  the  St.  John's  Sand  excluding  test  no. 

76  is  1104,  which  nearly  equals  the  value  for  the  dry  sand.  From  this  it  was 

concluded  that  for  sands  with  less  than  15  percent  fines  R  =  1100  could  be  used 

for  all  percents  saturation,  as  shown  by  the  first  of  equation  3. 

The  value  of  R  defines  the  relationship  between  TV  „  and  G  .  As 

r  max  max 

shown  in  Appendix  I,  equation  11  can  be  rearranged  to  eliminate  Tm  for 

reference  strain,  and  obtain  the  reference  strain  in  terms  of  G  as  given 

max 

by  equation  1.  Therefore,  the  data  in  figures  13  and  14  and  in  table  5  defining 
the  value  of  R,  are  the  basis  for  the  equations  given  in  Section  II,  paragraph  2 
for  determination  of  the  reference  strain,  and  indicate  the  degree  of  accuracy 
that  may  be  expected  when  using  that  part  of  the  practical  procedure. 
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Test 

No. 


Value 
of  R 


34 

1045 

35 

1094 

36 

1046 

37 

1026 

38 

1101 

39 

1116 

83 

990 

89 

1167 

101 

1125 

112 

loll 

Avg. 

1102 

3.  Comparison  of  Shear  Modulus  Values 

The  variations  of  normalized  shear  modulus  with  normalized  strain,  as 
given  by  equations  8  and  9,  are  shown  in  figure  15  for  different  values  of  a.  When 
a  =  0  the  relationship  is  hyperbolic.  Figure  15  shows  that  higher  values  of  a 
give  lower  values  of  G/G max  for  a  given  y/yr.  Using  measured  values  for  G, 
Gmax,  Tmax’  anc*  experimental  values  of  the  normalized  shear  modulus  and 
normalized  strain  can  be  calculated.  In  figures  16  through  22  such  experimental 
values  are  compared  to  the  values  given  by  equations  8,9,  and  10  of  the  practical 
procedure.  Figures  16  through  22  are  for  the  first  cycle  of  loading.  One  test 
of  each  different  soil  was  chosen  to  show  the  comparison.  The  tests  shown  were 
not  chosen  because  they  gave  the  best  comparisons.  The  parameters  necessary 
for  obtaining  the  calculated  curve  from  equations  8,9,  and  10  are  shown  in  the 
upper  right  hand  corner  of  each  graph.  The  shapes  of  the  calculated  curves 
correspond  to  various  values  of  a  between  0. 33  and  5. 34.  The  accuracy  of  the 
practical  procedure  is  quite  good  considering  the  range  of  soil  types  and 
conditions  to  which  they  are  applied. 

Comparisons  for  the  10th  and  100th  cycles  are  shown  in  figures  23  and 
24,  respectively.  For  these  comparisons  the  measured  value  of  normalized 
strain  along  with  the  value  of  a  calculated  from  equation  10  were  used  to  cal¬ 
culate  yh  from  equation  9.  This  value  of  y^  was  plotted  versus  the  measured 
normalized  shear  modulus,  thus  locating  the  different  symbols  in  figures  23  and 
24  representing  the  measured  relationship.  The  calculated  relationship  shown 
by  the  curves  in  figures  23  and  24  is  given  by  equation  8  and  is  hyperbolic, 
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Figure  15.  Variation  of  normalized  shear  modulus  with  normalized 
strain  as  given  by  the  practical  procedure  for  various 
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Figure  16.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  first  cycle,  WES  Sand,  St.  John's  Sand,  and 
Air  Force  Silty  Sand. 
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figure  17.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  first  cycle,  Air  Force  Silty  Clay,  Vicksburg 
Loess,  and  Vane eburg. 
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Figure  18.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  first  cycle,  Allen,  Kentucky  55,  and  Longhorn. 
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Figure  19.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  first  cycle,  West  Virginia  Shale,  Virginia 
Clay,  and  Dover. 
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Figure  21.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  first  cycle,  San  Francisco  Clay,  Ellsworth, 
and  Cheeks. 
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Figure  22.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  first  cycle,  Nevada  Clay. 


since  the  hyperbolic  strain  was  used  as  abscissa.  The  comparison  between 
measured  and  calculated  values  for  the  10th  and  100th  cycles  is  Rood. 

The  hyperbolic  strain  is  used  in  fibres  23  and  24  so  that  data  for 
different  soils  with  different  values  of  S  and  T  can  all  be  shown  on  the  same 
graph.  When  G/G  is  plotted  versus  y/yr,  the  expected  relationship 
depends  on  the  values  of  N,  S,  T  and  soil  type.  A  different  calculated  curve 
would  be  needed  for  each  test,  as  was  done  in  figures  16  through  22  for  N  1 
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Figure  23.  Comparison  of  measured  and  calculated  values  of  shear 
modulus,  10th  cycle. 
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Comparison  of  measured  and  calculated  values  of  shear 
modulus,  100th  cycle. 


SECTION  rv 


MIXED  LOADING  AMPLITUDES  AND  REST  PERIODS 

1.  Objective 

For  the  majority  of  tests  discussed  herein  a  cyclic  shear  loading  was 
applied  continuously,  where  the  next  cycle  began  immediately  as  the  preceding 
cycle  ended,  and  the  loading  amplitude  was  constant  for  ail  cycles  as  was  the 
rate  of  loading  and  unloading.  Because  actual  air  traffic  involves  varying  length 
rest  periods  between  load  applications  and  a  mixture  of  light  and  heavy  loadings, 
it  was  desirable  to  assess  these  effects  by  conducting  tests  where  the  loading 
could  be  programmed.  Then  a  soil  specimen  could  be  subjected  to  rest  periods, 
with  zero  load,  between  cycles  of  loading  and  the  amplitude  of  loading  could  be 
different  for  successive  cycles  of  loading. 

2.  Loading  Programs  and  Recorded  Stress-Strain  Relations. 

Three  primary  loading  sequences  were  used  for  mixed  amplitude  studies. 
The  first  was  an  increasing  load  sequence  as  shown  in  the  top  of  figure  25.  In 
this  figure  the  ordinate  is  proportional  to  the  applied  shear  stress  and  the 
abscissa  is  time.  The  slope  of  the  lines  in  this  figure  are  proportional  to  the 
rate  of  loading.  The  loading  program  in  the  top  of  figure  25  produces  a  sequence 
of  loads  with  monotonically  increasing  amplitudes  but  with  a  constant  rate  of 
loading  and  unloading.  The  second  loading  sequence  shown  in  the  bottom  of 
figure  25  produces  a  large  first  load  with  subsequent  loads  decreasing  in  ampli¬ 
tude  monotonically.  The  third  sequence  was  a  combination  of  these  two  where 
i.ie  sample  was  subjected  first  to  the  increasing  load  sequence  followed  by  the 
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decreasing  load  sequence. 

Examples  of  the  recorded  stress-strain  relations  for  the  increasing 
load  sequence  are  shown  in  figures  26,  27,  and  28  for  silty  clay,  silty  sand, 
and  WES  loess  specimens,  respectively.  Figure  26a  shows  the  results  of 
applying  46  cycles  of  increasing  load  to  a  specimen  of  silty  clay.  These  46 
cycles  were  applied  continuously  without  a  rest  period.  At  the  end  of  the  46th 
cycle  there  was  a  short  rest  period  while  the  recording  paper  was  changed. 

After  the  rest  period  the  increasing  load  sequence  was  continued  and  the  results 
of  cycles  47  through  67  were  recorded  in  figure  26b.  The  recording  was 
approaching  the  top  of  the  paper  on  cycle  63  necessitating  a  change  in  stress  scale 
during  this  cycle.  The  scales  corresponding  to  various  sections  of  the  recording 
are  shown  on  the  figure.  At  the  end  of  cycle  67  there  was  a  short  rest  perijd 
as  the  recording  paper  was  changed.  The  increasing  load  sequence  was  again 
continued  and  the  results  for  cycles  68  through  83  are  shown  in  figure  26c.  At 
the  end  of  the  82nd  cycle  the  pen  of  the  recorder  was  shifted  electronically  to  the 
left  and  up  because  the  strain  was  becoming  large  with  the  recording  approaching 
the  edge  of  the  paper.  The  loading  curve  for  cycle  83  is  shown  above  the  main 
recording  and  as  can  be  seen  failure  of  the  specimen  occurred  on  the  83rd  cycle. 
Similar  procedures  were  used  for  the  recordings  in  figures  27  and  28. 

The  result  of  applying  the  decreasing  load  sequence  shown  in  the  bottom 
of  figure  25  to  a  specimen  of  silty  sand  is  shown  in  the  top  of  figure  29.  This 
same  decreasing  load  sequence  was  applied  to  a  specimen  of  silty  clay  producing 
the  results  shown  in  the  bottom  of  figure  29.  However,  the  specimen  of  silty 
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Figure  26.  Recorded  stress-strain  relation,  increasing  load  sequence 
silty  clay. 
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Recorded  stress-strain  relation,  increasing  load  sequence,  silty  sand 
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Figure  28,  Recorded  stress-strain  relation,  increasing  load  sequence 
Vicksburg  Loess. 
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Figure  29.  Recorded  stress-strain  relation,  decreasing  load  sequence. 
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clay  had  previously  been  subjected  to  the  increasing  load  sequence  in  the  top  of 
figure  25.  From  these  results  it  can  be  seen  that  for  decreasing  load  sequences 
there  is  little  additional  accumulation  of  set  strain  as  the  smaller  and  smaller 
loads  are  applied.  Hence,  the  curves  for  successive  loads  trace  over  the  same 
area  of  the  recording. 

3.  Effect  on  the  Shear  Modulus 

The  values  of  normalized  modulus  for  all  of  the  mixed  amplitude  tests 
are  plotted  versus  normalized  strain  for  the  silty  clay  in  figure  30  and  for  the 
silty  sand  in  figure  31.  The  open  symbols  are  for  increasing  amplitude  loading- 
sequences  and  the  solid  symbols  are  for  decreasing  amplitude  loading  sequences. 
Also  shown  on  these  figures  are  the  curves  calculated  from  equations  8  and  9  for 
values  of  a  =  0  and  a  =  2.  As  shown  in  figure  17,  the  value  of  a  for  the  first 
cycle  of  test  no.  21  on  the  Air  Force  silty  clay  is  3. 50.  This  value  decreased 
with  number  of  cycles  and  approaches  zero  for  a  large  number  of  cycles.  For 
the  first  cycle  of  test  no.  4  on  the  Air  Force  silty  sand  in  figure  16  the  value  of 
a  =  5.34.  Again  this  value  approaches  zero  for  a  large  number  of  cycles.  Since 
most  of  the  values  shown  in  figures  30  and  31  fall  in  the  range  between  a  =  0 
and  a  =  2 ,  it  appears  that  the  effects  of  the  mixed  amplitudes  and  rest  periods 
on  the  shear  modulus  are  smaller  than  the  scatter  in  measured  values  from  test 
to  test  of  the  same  soil.  Therefore,  it  is  concluded  that  for  practical  purposes 
the  procedure  outlined  in  Section  II  can  be  applied  to  mixed  traffic  conditions. 
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Figure  31.  Normalized  shear  modulus  versus  normalized  strain  for 
silty  sand,  mixed  amplitudes  and  rest  periods. 


SECTION  V 


CONC  LUSIONS 

1.  A  practical  procedure  for  reducing  the  shear  modulus  of  soils  with 
increasing  strain  amplitude  has  been  developed.  It  has  been  shown  that  the 
procedure  gives  reasonably  accurate  results  compared  to  values  measured  in 
the  laboratory,  for  a  wide  variety  of  soil  types  and  conditions.  The  study  of 
mixed  amplitudes  and  rest  periods  indicates  that  th~*  procedure  can  also  be  applied 
to  mixed  traffic  conditions. 

2.  The  shear  modulus  of  soils  in  general  varies  between  extreme  limits. 

It  may  commonly  be  less  than  1000  psi  or  greater  than  50000  psi.  Likewise,  the 

relationship  between  shear  modulus  and  strain  level  is  extremely  variable.  The 

picture  can  be  simplified  by  normalization.  The  shear  modulus  is  normalized  by 

considering  G/G  ,  but  the  relationship  between  normalized  shear  modulus 
max 

and  strain  level  is  still  quite  variable;  because,  a  given  sti-ain  does  not  have  the 
same  effect  on  all  soils  or  on  the  same  soil  under  different  states  of  stress. 

When  strain  is  also  normalized  by  dividing  by  the  reference  strain,  yr  Tmax/ 
Gmax»  a  single  relationship  for  all  soils  between  normalized  shear  modulus 
and  normalized  strain  is  approached.  However,  this  relationship  is  still 
affected  somewhat  by  soil  type,  percent  saturation,  and  number  of  cycles  and 
rate  of  loading.  A  single  relationship  for  all  soils  and  conditions  is  finally 
obtained  by  defining  the  hyperbolic  strain  that  depends  mainly  on  the  normalized 
strain  but  accounts  for  the  residual  effects  of  soil  type,  S,  N,  and  T.  Hence, 
these  concepts  of  reference  strain  and  hyperbolic  strain  have  made  possible  the 
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development  of  the  practical  procedure  of  paragraph  1. 

3.  The  123  simple  shear  tests  of  24  different  soils  reported  herein  show  the 
following  parameter  effects  on  the  relationship  between  normalized  shear  modulus 
and  normalized  strain.  For  a  given  value  of  normalized  strain,  the  normalized 
shear  modulus  increased  with  number  of  cycles  of  loading,  decreased  with 
increasing  percent  saturation,  and  increases  with  rate  of  loading.  There  appears 
to  be  very  little  effect  of  void  ratio  on  this  relationship. 

4.  The  reference  strain  that  is  so  important  to  the  procedure  of  paragraph 
1  depends  on  r  and  G  .  For  pavement  evaluation  the  value  of  G  is  to 
be  determined  by  the  nondestructive  vibratory  test.  It  is  not  desirable  to  have  to 
measure  the  value  of  Tmax<  In  this  study  a  relationship  between  reference  strain 
and  G  has  been  established,  eliminating  the  necessity  for  estimating  r 

mux  iiitu/v. 

The  relationship  depends  on  void  ratio,  percent  saturation  and  plasticity  index 
of  the  soil.  It  is  best  defined  for  nonplastic  and  low  plasticity  soils,  where  the 
data  in  figures  13  and  14  indicate  that  the  error  is  less  than  25  percent  for 
about  50  percent  of  the  cases.  It  should  be  remembered  that  a  25  percent  error 
in  the  value  of  yr  does  not  produce  a  correspondingly  large  error  in  the  deter¬ 
mination  of  G.  For  a  value  of  normalized  strain  of  one,  a  25  percent  error  in 
reference  strain  produces  about  12  percent  error  in  the  determination  of  G. 

The  relationship  is  much  less  well  defined  for  high  plasticity  soils  with  liquid 
limit  greater  than  50,  where  the  error  in  determination  of  reference  strain  may 
be  greater  than  100  percent. 

5.  Assuming  the  normalized  strain  is  accurately  known,  the  error  to  be 
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expected  in  the  relationship  between  normalized  strain  and  normalized  shear 
modulus  is  shown  by  the  data  in  figures  16  through  24,  where  the  vast  majority 
of  points  fall  within  plus  or  minus  20  percent. 

G.  Finally,  an  electromagnetic,  hollow  cylinder,  simple  shear  test  capable 
of  accurate  measurement  of  the  stress-strain  relations  for  soils  over  a  wide 
range  of  strains,  from  about  10“5  in/in  to  strains  corresponding  to  failure,  and 
where  the  loading  can  be  programmed  to  test  effects  of  mixed  amplitudes  and 
rest  periods  has  been  developed. 


APPENDIX  I 


DERIVATION  OF  REFERENCE  STRAIN  EQUATION 
Neglecting  the  effects  of  overconsolidation,  for  a  simple  shear  state  of 


stress 


Tmax  ~  ao  sin  0  (12) 

where  5a  =  effective  mean  principal  stress  and  0  =  effective  angle  of  shearing 
resistance.  Also  as  stated  in  reference  1, 

-  1/2 

(13) 


-  1/2 

G  =  1230  F  or 
max  ° 


Equation  13  is  based  on  many  tests  of  saturated  soils.  Solving  for  aQ  in 
equation  13,  substituting  into  equation  12  and  re.  lacing  the  factor  1230  by  It, 


since  this  factor  may  not  be  constant  for  partially  saturated  soils 


^max  sin  * 


max 


(14) 


As  shown  in  reference  2 ,  the  effective  angle  of  shearing  resistance  can  be 
defined  approximately  by 


sin  0  "0.6  ~  0.25  (PI) 


0.0 


(15) 


Substituting  equation  15  into  equation  14 


T  =£jnax  [0.6  -  0.25  (PI)0'6] 
max  r2  p2 


(16) 


The  reference  strain 


y  = 

r  g 


max 


(17) 


max 


2. 


Lambe,  T.  W.  and  Whitman,  It.  V. ,  Soil  Mechanics,  Wiley,  1969, 
p.  307. 
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Substituting  equation  16  into  equation  17 

y  [0.6  -  0.25  (PI)  ° 

r  RY 
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APPENDIX  II 


TESTING  METHODS  AND  PROCEDURES 

A  simple  shear  test,  where  a  hollow  cylinder  of  soil  was  confined  »n  a 
pressure  chamber  and  loaded  torsionally  about  the  axis  of  the  cylinder,  was 
employed  for  this  research.  The  torsional  load  was  applied  electromagnetically 
by  reguiating  the  voltage  to  four  large  coils  within  the  fields  of  four  corresponding 
permanent  magnets.  The  system  was  capable  of  applying  a  maximum  torque  of 
approximately  60  kg-cm.  The  voltage  was  produced  by  a  wide  frequency  range, 
triangular  wave  signal  generator,  amplified  by  a  50  watt  DC  amplifier.  With 
this  electromagnetic  system  the  rate  and  amplitude  of  loading  could  be  changed 
by  simply  turning  the  appropriate  knob  on  the  signal  generator.  The  rate  of 
loading  for  the  tests  varied  from  approximately  0.2  to  450  kg/cm  per  hr. 

Inertial  forces  were  negligible  for  these  rates  of  loading.  Torque  and  angular 
motion  were  measured  with  electrical  transducers  located  inside  the  pressure 
chamber,  thus  eliminating  the  influence  of  friction  and  apparatus  deformation  on 
the  measurements.  The  signals  from  these  transducers  were  recorded  with  an 
X-Y  recorder. 

Most  of  the  soil  specimens  tested  were  subjected  to  approximately  1000 

cycles  of  a  constant-amplitude  shear  stress.  Examples  of  the  recorded  shear 

stress-strain  relations  are  shown  in  figure  32.  This  figure  illustrates  the  range 

of  strain  amplitudes  tested.  For  test  no.  7,  top  of  figure  32,  the  cyclic  shear 

stress  was  approximately  80  percent  of  the  strength  of  the  sample.  For  test 

no.  15,  record  no.  1,  bottom  of  figure  32,  the  cyclic  shear  stress  was  only  2.5 
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0.00062  % 
0.0021  kg/cm^-f 


SILTY  SAND 
e=  0.45 
S  =  100°/o 

^ham =  0 k9/cm 
9  sec /cycle 


Figure  32.  Typical  recorded  stress-strain  relations 


percent  of  the  strength  of  the  sample.  The  entire  graph  sheet  for  test  no.  15,  if 

changed  to  the  scale  of  test  no.  7,  would  fit  in  the  small  black  rectangle  at  A  in 

the  lower  left  corner  of  test  no.  7.  Cyclic  loading  tests  were  conducted  for 

strain  amplitudes  as  small  as  about  10~5  in/in  to  strain  amplitudes  as  large  as 
-2 

0. 5  x  10  in/in.  After  the  cyclic  loading  of  each  sample,  the  load  was  increased 
to  failure  in  order  to  measure  T 

max 


63/64 


